A multivariate multifractal model for return fluctuations 
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4L i Abstract 

c^ . 
-)— > , 

j^ I In this paper we briefly review the recently inrtroduced Multifractal Random Walk (MRW) that 

^ ' is able to reproduce most of recent empirical findings concerning financial time-series : no correlation 

' . between price variations, long-range volatility correlations and multifractal statistics. We then focus on 

^ I its extension to a multivariate context in order to model portfolio behavior. Empirical estimations on real 

O ' data suggest that this approach can be pertinent to account for the nature of both linear and non-linear 

1 ^i ' correlation between stock returns at all time scales. 

T-H ! 

> ■ 1 Introduction 

O; 

^vq I Multifractal processes and the deeply connected mathematics of large deviations and multiplicative cas- 

0^ ■ cades have been widely used in many contexts to account for the time scale dependence of the statistical 

^r , properties of a time-series. Recent empirical findings ||l|, Q S H^ suggest that in finance, this framework 

f"^ I is likely to be pertinent. The recently introduced Multifractal Random Walks (MRW) [^ are multifractal 
processes that have proved successful to model return fluctuations. They can be seen as simple "stochas- 
tic volatility" models (with stationary increments) whose statistical properties can be precisely controlled 

C I across the time scales using very few parameters. In that respect, they reproduce many features that 

'^ ■ characterize market price changes [|l^] including the decorrelation of the price increments, the long-range 

G . correlation of the volatility and the way the probability density function (pdf) of the price increments 

fi I changes across time-scales, going from quasi Gaussian distributions at rather large time scales to fat tail 

L" ' distributions at fine scales. 
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In a recent work |11], Muzy et. al. have elaborated a "multivariate multifractal" framework that 
accounts for the time scale dependence of the mutual statistical properties of several time-series. They 
C^ '_ have shown that the statistical properties of financial time-series can be described within this framework. 
Though initially introduced for modelling single asset variations, the MRW models can be naturally 
extended in order to fit this new multivariate framework. Thus, the so-obtained Multivariate MRW 
(MMRW) can be used to reproduce precisely the statistical properties of several assets at any time-scale. 
This is of course particularly useful for modelling a portfolio behavior. 

The goal of this paper is to explain how MMRW models are built and to show, using real data, 
that eventhough they involve very few parameters, they allow one to capture not only linear correlation 
between assets but also non linear correlation at all time scales. The paper is organized as follows. In 
section 2, after recalling briefly the main notations and definitions involved in the "classical" monovariate 
multifractal framework, we introduce the MRW model defined in [g], recall its main properties and 



show that it reproduces precisely, at any time-scale, the statistical properties of real financial time- 



series. In section 3, after presenting the multivariate multifractal framework defined in |11|, we introduce 
the MMRW model and perform some analytical computations of its multifractal statistics. Numerical 
estimations of the parameters in the case of financial times-series are extensively discussed. Conclusions 
and prospects are reported in section 4. 

2 The Multifractal Random Walk (MRW) model 

2.1 Multifractal processes and cascade models 

A multifractal process is a process wich has some scale invariance properties. These properties are 
generally characterized by the exponents Qq which govern the power law scaling of the absolute moments 
of its fluctuations, i.e., 

M{q,l)=Kqf^, (1) 

where 

M{q, l)=E {\5iX{t)\'i) = E {\X{t + 1)- X(i)r) , 

where X{t) is supposed to be a stochastisc process with stationary increments. Some very popular 
stochastic processes are the so-called self-similar processes [|l5|. They are defined as processes X{t) 
which have stationary increments and which verify (in law) 

5iX{t) ='""' {l/L)"5LX{t), V/,L > 0. 

For these processes, one easily gets Cq = qH^ i.e., the C,q spectrum is a linear function of q. Widely used 
examples of such processes are (fractional) Brownian motions (fBm) or Levy walks. 

However, many empirical studies have shown that the C,q spectrum of return fluctuations is a non 
linear convex function. Let us note that, using a simple argument, it is easy to show that if C,q is a non- 
linear convex function the scaling behavior (||) cannot hold for all scales I but only for scales smaller than 
an arbitrary large scale T that is generally referred to as the integral scale. A very common approach 
originally proposed by several authors in the field of fully developed turbulence |12, 14, P, ^, |5|, has 



been to describe such processes in the scale domain, describing the cascading process that rules how 
the fluctuations evolves when going from coarse to fine scales. Basically, it amounts in stating that the 
fluctuations at the integral scale T are linked to the ones at a smaller scale / < T using the cascading 
rule 

5iX{t) ='»- Wi,T5TX{t) (2) 

where Wu^ is a log infinitely divisible stochastic variable which depends only on the ratio l/T. A 
straightforward computation pi then shows that the pdf Pi{5X) of 5iX changes when varying the time- 
scale / according to the rule 

Pi{5X) = I Gi/T{u)e-^PT{e-^dX)du, (3) 

where the self-similarity kernel Giij- is the pdf of In Wnx- Since Wnx is a log infinitely divisible variable, 
the Fourier transform of Gi/t is of the form 

Gi/Tik) = &^'/'^ik). (4) 

From that equation, one easily gets the expression of the Cq spectrum 

Cq= In G{-iq). (5) 

Thus, the simplest non-linear case is the so-called log-normal model that corresponds to a parabolic Cq 
and a Gaussian kernel. 



Multiplicative cascading processes are examples of processes satisfying the cascading rule (^). How- 
ever, they have fundamental drawbacks: they do not lead to stationary increments and they do not 
have continuous scale invariance properties, i.e., Eq. (^ and consequently Eq (|^) only holds for discrete 
scales /„ = A". To our knowledge, the MRW's are the only known multifractal processes with continuous 
dilation invariance properties and stationary increments. 

2.2 Introducing the MRW model 

An MRW process X{t) is the limit process (when the time discretization step At goes to 0) of a standard 
random walk XAi[A;] with a stochastic variance (volatility), i.e., 

X(t) = lmiXAi(t), 

with 

t/At 

XAt{t) = Y. eAt[k]e^^'^''\ 

k=l 

where e'^'^'^^''' is the stochastic volatility and e^t a gaussian white noise of variance a'^At and which is 
independant of u!/\f The choice for the process wa* is simply dictated by the fact that we want the 
scaling (IT]) to be exact for all time scales / < T. Some long but straightforward computations show 
that this is achieved if wa* is a stationary Gaussian process such that -©(wAti^]) = — Var (lja*!^]) a-nd 
whose covariance is 

Cov{u;At[k],uJAt[l]) = X'^\npAt[\k - l\] 

where 

f n^tL\A, for \k\ < T/At - 1 

PAt[k] = { fi+i)^* .' 'r ' 

I 1 otherwise 

Let us note that it corresponds to a log-normal volatility which is correlated up to a time lag T. 
One can then prove |Q] the multifractal scaling property 

M{q,l)=Kgf\ yi<T, (6) 

with 

Q = {q-q{q-2)X^)/2. (7) 

Since Cq is a parabolic function, it indicates that the self-similarity kernel Gwj^ which links the pdf's at 
different time scales (Eq. (|^) is Gaussian. Moreover one can show |^] that the magnitude correlation 
Cui{1,t) defined by 

C^(At, = Cov (In \SiX{t)\,ln \6iX{t + r)|) , (8) 

behaves like 

9. /At 



C^(At,/)~-A^ln( — ), KT. (9) 



2.3 Modelling return fluctuations using MRW 



MRW processes can be used to model return fluctuations ||l^. For this purpose 3 parameters need to be 
estimated : the variance a, the integral scale T and the intermittency parameter A (it is called this way 
since it controls the non linearity of the C,q spectrum and consequently it controls "how much stochastic" 
is the variable Wux)- The variance a can be estimated using the simple relation Var{X{t)) = a'^t. Both, 
the decorrelation scale T and the parameter A can be obtained from the expression (^ of the magnitude 
correlation. Let us note that A can be also estimated independantly from the (g spectrum (Eq. (0)). The 
consistency between these two completly different estimators of A is a very good test for the validity of 
the model. 
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Figure 1: Modelling intraday Japenese Yen futures using MRW. The MRW is used to model 
the de-seasonalized logarithm of the time-series displayed in (a). The parameters have been estimated to 
a^ = 4.10"^, A^ = 0.03 and T = 1 year, (a) Plot of the original index time-series : Japenese Yen futures 
from March 77 to February 99 (intraday tick by tick data), (b) Plot of a sample time series of length 
2^^ of the model, (c) Cg spectrum estimations for the Yen futures fluctuations (o) and for the MRW 
model (x). The solid line corresponds to the theoretical prediction (Eq. (|^). (d) Magnitude correlation 
function estimations as defined in Eq. (^) (with / ~ 4 days) for the Yen futures fluctuations (o) and for 
the MRW model (x). The solid line corresponds to the theoretical prediction (Eq. (P)). 



Parameter estimations have been made on financial data (japenese Yen futures). As shown in figure 
1, the MRW reproduces very precisely both the parabolic Cg spectrum (which describes, through Eqs 
®5 (§)' ^^^ ^^^ return fluctuation pdf evolves when going from one time scale to another) and the 
correlation structure of the magnitude. 

Let us remark that one can show that Kg = -|-oo (in Eq. ^) if (^g < 1 and thus the pdf of 6iX{t) has 
fat tails Q. In order to control the order of the flrst divergent moment (without changing A), one could 
simply choose for the eAt's a law with fat tails (e.g. t-student laws). 



3 The multivariate multifractal random walk (MMRW) model 

3.1 The multivariate multifractal framework 

In this section, we generalize in a very natural way the multifractal framework introduced in section 2.1 
to multivariate processes. This generalization is inspired from ||Il|. It basically consists in rewriting the 
cascading rule (|2|) using multivariate processes. Thus if X = {Xi, . . . ,Xj\f) is a multivariate process, we 
will assume that it satisfies 



{SiXiit}} 



_law 



l<i<N 



{W,^l,T5TX,{t)] 



l<i<N 



where W = < Wi^i/x r . is a log infinitely divisible stochastic vector which depends only on the ratio 

l/T. A straightforward computation [^ then shows that the pdf Pi{5'K) of 5/X changes when varying 
the time-scale / according to the rule 

Pi{6X) = f dui... /"(iuivG,/r(u)e-"i---'^^Pr(e""0(5X), (10) 

where u = {ui, . . . ,uj^} and where the Gi/t is the pdf of lnW;/y. (We used the notation a (g) b = 
(ai6i, . . . , OAr6jv)). In the same way as in section 2.1, one can easily get the scaling law of the moments 

M(gi, ...,qN,l)=E {\6iXi{t)\'^' . . . {SiX^it)]'^^) = Kg.^gJ^"^ "n , (11) 

where the multifractal spectrum Cqi,...,qN ^^ linked to the self-similarity kernel through the relation 

Cqi,...,qN = G{-iqi,...,-iqN), 
*equation where G is defined as in Eq. (Q). 

3.2 Introducing the MMRW model 

In order to account for the fluctuations of flnancial portfolios and to consider management applications 
of our approach, it is important to build a multivariate version of the MRW model. Since only Gaussian 
random variables are involved in the construction of section 2.2, this generalization can be done in a very 
natural way 0. The MMRW walk X(t) is defined as 

t/At 
X(t) = limXAi(t),= lim V €At[k] e'^^'^. 

k=l 

(We again used the notation a ® b = (ai6i, . . . ,a7v^Ar))- The process ea* is Gaussian with zero mean 
and covariance Cov{ei^At{t), ej^^tit + t)) = (5(r)I]jjAi. The matrix S quantifies the variance and the 
correlation of the different white noises involved in each component of X. We will refer to this matrix as 
the "Markowitz matrix". The magnitude process to At is Gaussian with covariance Cov{uJi^At{t),^j,AT) = 
Aij ln(Tij / {At + |r|)) (for At + \t\ < Tij) and elsewhere, where the matrix A controls the non-linearity 
of the multifractal spectrum so we will refer to it as the "multifractal matrix". Moreover, as in section 
2.2, the mean of the process is chosen so that -^(wAt) = —Var(iJ^t). Let us note that the previously 
defined coefficients o"^ and A^ for an asset i correspond respectively to the diagonal elements Sjj and An. 
In order to show that MMRW are multivariate multifractal processes (within the framework of the 
previous section), we would like now to compute the C(ji,...,g]v spectrum. There are 2 cases for which 
this computation is basically the same as for the regular MRW model : (i) the case were all the white 
noises are decorrelated (i.e., Xl is diagonal), (ii) the case where the stochastic variances of all the assets 
correspond to the same process, i.e., LUi^At = ^jAti "^hj- In both cases, a straightforward computation 
shows that the scaling law ( [ll| ) holds VZ < minjj(Tjj) and the spectrum is 

N 

Cqu...,qN=Y.Q- Yl ^ijQiQj: (12) 

i=l l<i<j<N 

where C refers to the spectrum of the Xi component of X. The computation of the spectrum is trickier in 
the general case. However one can show that, in this case, all the extra terms (compared to the particular 
case (i)) that appear in the development of M{qi, . . . , q^, I) go to when At ^ 1^. Consequently, the 
multifractal spectrum has the same expression. 

Since the spectrum is a parabolic function, it indicates that the self-similarity kernel Gi/x which links 
the pdf's at different time scales (Eq. (|^)) is Gaussian. Moreover one can show [^ that the magnitude 
correlation behaves like 

Gov (In \6rXi{t)\,ln \6rXj{t + /)|) ~ -Aij ln(/) + C . (13) 
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Figure 2: Estimation on real data of a non diagonal term Ajj of the multifractal matrix A. 

Ai2 is estimated using daily data (from 1993 to 2000) where the first asset is AXA and the second one is 
TOTAL-FINA. As explained in the text, two estimators can be used. The solid curve corresponds to the 
Ri2,q=i estimator (Eq. (15)) which leads to the estimation A12 ^ 0.032. The dashed curve corresponds 



to the covariance estimator (Eq. (O)) which leads to the estimation A12 — 0.036. These two estimators 
lead to consistent estimation of A12. 



3.3 Modelling return fluctuations using MMRW 

For modelling a basket of assets using an MMRW, one needs to estimate the Markowitz matrix S, the 
multifractal matrix A and the different integral scales Tij. As for the MRW, the magnitude correlation 
can be used for estimating both the integral scales and the multifractal matrix. The Markowitz matrix 
S can be estimated, after having estimated A, by using the simple relation 



Cov{Xiil),X,il)) 



MJ^u- 



Uj)l ^ 



(14) 



Let us note that A can be also estimated independantly from the C(ji,.-,g]v spectrum (Eq. (|^)). Indeed, 
to estimate Ajj one could simply estimate, for instance, the exponent of the power law scaling 



*i)QVv 



E{\X,m)E{\XM') 



l-A,jq'^ 



(15) 



The consistency between these two different estimators of A is a very good test for the validity of the 
model. 

Parameter estimations have been made on some assets (daily data) of the cac40. Figure 2 shows 
the estimations of the non diagonal term A12 of the multifractal matrix A using both the magnitude 
correlation estimator (Eq. (|13|)) and the i?i2,g(Af) estimator (Eq. (|l^)). These two different estimators 
lead to very close estimations of A12 which is consistent with the MMRW model. 

On Figure 3 we have displayed the histograms of all the non diagonal terms Ajj (resp. Sjj and Tij) 
for all the pair of assets in the cac40. Eventhough the histogram of the Aij is pretty wide, its maximum 
is reached for Aj, ~ 0.02 which is the most common value found when estimating An = A? in the 



monovariate case |1C]. In the same way the histogram of the Tij (resp. I^jj) has a peak around 1-2 years 
(resp. 4.10"^) which also corresponds to the most common value found when estimating Tj (resp. a^) 
in the monovariate case [I^. These results suggest that, as a first approximation, one could model these 
assets using an MMRW which shares the same magnitude process for all the assets, i.e., tOi = ujj V i,j. 
Though it is clearly not exactly the case, it simplifies the model a lot and allows to perform many analytic 
computations. 
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Figure 3: Parameters estimations for a basket made of all the cac40 assets. Histogram of the 
estimations, for all the pair of assets in the cac40, of the non diagonal terms (a) Ajj, (b) Tij and (c) "Sij 
(using Eq. |^). 



4 Conclusion 

In this paper, we have introduced a multivariate model for return fluctuations that is a definite step 
beyond standard correlation analysis. It basically corresponds to a multivariate random walk using a 
stochastic volatility. This model can be characterized using the recently introduced notion of multifractal 
multivariate that itself relies on the idea that the simplest way to describe the statistics of a process at all 
time scales is to assume some scale invariance properties. Consequently, the MMRW has a potentiality to 
capture the whole return joint law of a basket of assets at all time horizons. As shown in this paper, we 
are able to reproduce the main observed characteristics of financial time-series: no correlation between 
price variations, long-range volatility correlations, linear and non-linear correlation between assets and 
the price increment pdf and the way it changes when varying the time-scale. All of these features 
can be controlled using a few parameters : the multifractal matrix which controls the scale invrariance 
properties, the integral scales which controls the volatility correlation and the Markowitz matrix which 
controls the noise correlation. Moreover, as we have already pointed out, in good approximation, one 
can consider that all the assets share the same volatility process. Not only it reduces the number of 
parameters but it makes any analytical computation much easier. We are currently applying MMRW for 
portfolio management and risk control. 
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